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Abstract 

In this paper we study the pathwise stochastic Taylor expansion, in the sense of our previous 
work [3] , for a class of Ito-type random fields in which the diffusion part is allowed to contain 
both the random field itself and its spatial derivatives. Random fields of such an "self-exciting" 
type particularly contains the fully nonlinear stochastic PDEs of curvature driven diffusion, as 
well as certain stochastic Hamilton- Jacobi-Bellman equations. We introduce the new notion 
of "ra-fold" derivatives of a random field, as a fundamental device to cope with the special 
self-exciting nature. Unlike our previous work [3], our new expansion can be defined around 
any random time-space point (r, £), where the temporal component r does not even have 
to be a stopping time. Moreover, the exceptional null set is independent of the choice of 
the random point (r, £). As an application, we show how this new form of pathwise Taylor 
expansion could lead to a different treatment of the stochastic characteristics for a class of 
fully nonlinear SPDEs whose diffusion term involves both the solution and its gradient, and 
hence lead to a definition of the stochastic viscosity solution for such SPDEs, which is new in 
the literature. 
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1 Introduction 

In our previous work [3] we studied the so-called pathwise stochastic Taylor expansion for a class 
of Ito-type random fields. The main result can be briefly described as follows. Suppose that 
{u(t,x), (t,x) 6 [0, T] x MJ 1 } is an Ito-type random field of the form 



where B is a 1-dimensional standard Brownian motion, defined on a complete probability space 
(ft, T , P). If we denote F = {J^ } t >o to be the natural filtration generated by B and augmented by 
all P-null sets in T , then under reasonable regularity assumptions on the integrands u\ and U2, the 
following stochastic "Taylor expansion" holds: For any stopping time r and any ^-measurable, 
square-integrable random variable £, and for any sequence of random variables {(t^,^)} where 
Tfc 's are stopping times such that either > t, t^t; or < t, t^t, and 's are all T Tk i\ T - 
measurable, square integrable random variables, converging to £ in L 2 , it holds almost surely that 



where (a,b,c,p,q,X) are all T ' T -measurable random variables, and the remainder o(Cfc) are such 
that o(Ck)/Ck — > as k — > oo, in probability. Furthermore, the six-tuple (a,b,c,p,q,X) can be 
determined explicitly in terms of u\, u<i and their derivatives. 

By choosing u\ and U2 in different forms, we then extended the Taylor expansion to solutions 
of stochastic differential equations with initial state as parameters, and to solutions of nonlinear 
stochastic PDEs. In the latter case we further introduced the notion of the stochastic super- (sub- 
)jets using the Taylor expansion, from which the definition of stochastic viscosity solution was 
produced. We should note that in [3] all the SDEs and SPDEs have the diffusion coefficient in the 
form g(t, x, u(t, x)), that is, they only involve the solutions themselves. Such a structure turns out 
to be essential for the so-called Doss-Sussmann transformation, and in that case the stochastic 
viscosity solution became natural. 

In this paper we are interested in the stochastic Taylor expansion for random fields of the 
following form: 



u(t, x) = u(Q, x) + / m(s,x)ds+ (g(t,x,Du(s,x)),dB s ), (t, x) 6 [0, T] xl™, (1.3) 



where u\ is a random field, and g is a deterministic function. We shall assume that they are 






+ (|r fe -r|) + (|e fe -e| 2 ) 
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"smooth" in the sense that all the desired derivatives exist, and the degree of smoothness will be 
specified later. 

Given such a random field, we again consider the possibility to expand u in the sense of (jl.2p . 
but this time in a more natural way: We shall allow the pair (r, £) : $7 — > [0, T] x to be 
arbitrary random points. Furthermore, we should note that although the Taylor expansion (jl.2p 
holds almost surely, in general the null set may depend on the choice of random point (t, £). In 
this paper we shall look for a universal expansion, in the sense that there is a subset $7 C SI with 
full probability measure, on which the stochastic Taylor expansion holds for all choices of random 
points (r, £). These technical improvements make the stochastic Taylor expansions much more 
"user- friendly" , and more importantly, it will be more effectively used in our study of stochas- 
tic viscosity solution, especially in the proof of uniqueness, as we shall see in our forthcoming 
publications on that subject. 

As one can easily observe from (jl.3p that the random field u is actually already a solution 
to a first order stochastic PDE. An immediate difficulty in deriving the Taylor expansion is the 
characterization of the derivatives of the random field, as it will increase in a "bootstrap" way, 
very similar to those that one has often seen in the anticipating calculus. As a consequence, the 
original Doss-Sussmann type transformation used in our previous works [H [SJ E] no longer works 
in this case. To overcome this difficulty we introduce the notion of "n-fold derivative" of a random 
field, which essentially takes (u, Du, ■ ■ ■ , L> n_1 u) as a vector-valued random field, and define its 
derivative in a recursive way. Such a definition turns out to be very close to the idea of converting 
a higher order ordinary differential equation to a first order system, and is mainly motivated by 
the "stochastic characteristics" of a stochastic PDE (cf. e.g., Kunita [12J). In fact, by combining 
the definition of stochastic characteristics in [12J and the stochastic Taylor expansion developed in 
this paper, we are able to rigorously define a stochastic diffeomorphism that relates the Stochastic 
PDE of the form 

du(t,x) = f{x,u,D x u,D 2 x u)dt + g(x,D x u)) dBl (t, x) € (0, T) x R d ; 
u(Q,x) = u (x), xeK d , 
to a PDE without Brownian components in their Taylor expansions around any temporal-spatial 
point, generalizing the Doss-Sussmann transformation in our previous works to the present case. 

We should note that one of the main difficulties in the study of the stochastic viscosity solution 
can be described as "local" vs. "global". That is, the local nature of the viscosity solution vs. 
the global nature of the stochastic analysis (e.g., stochastic integrals). Our idea is to "localize" 
the stochastic integral, or solution to the SPDE via the stochastic Taylor expansion that we 
established in the previous section. It should be noted that the universal set "fi"' that we found 
in these expansions is the essential point here. 
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Finally, we would like to point out that the topic of stochastic Taylor expansion has been ex- 
plored by many authors in various forms, and used to provide numerical and other approximation 
schemes for SDEs and SPDEs or randomized ODEs and PDEs (see, for example, [T| [3"1 1101 to 
mention a few). These expansions often use either the Lie-algebraic structure of the path space 
or the chaotic type expansion of multiple stochastic integrals. As a consequence it is hard to de- 
duce the simultaneous spatial-temporal expansions that we are pursuing in this work, especially 
when the remainders are estimated in a pathwise manner. We should also note that the path- 
wise version of stochastic viscosity solutions, suggested by Lions-Souganidis [El E], has found 
an effective framework recently, using the theory of rough path (cf. Caruana-Friz-Oberhauser [7]). 
However, due to the special nature of the rough path integrals, the arguments seem to depend 
heavily on the fact that there exist stochastic characteristics in the form of C 3 -diffeomorphisms 
which transform the SPDE to a pathwise PDE. Consequently, the SPDE studied in [7], while fully 
nonlinear in the drift, seems to be restricted to the cases when the diffusion coefficient depends 
only (linearly) on Du, the gradient of the solution, so that a chain-rule type of argument could 
be applied. The generality of the diffusion part in the fully non-linear PDE suggested in this 
paper, and the stochastic characteristics related to it, does not seem to be an easy consequence of 
such a method. Furthermore, our Taylor expansions are constructed within a more "elementary" 
stochastic analysis framework, by exploiting the properties of Brownian motion without using the 
advanced algebraic geometric structure of the path spaces, therefore we believe that it provides a 
more accessible alternative. 

This paper is organized as follows. In Section[2]we clarify all the necessary notations, and state 
the main theorem. In section 3 we give a fundamental estimate of this paper, regarding multiple 
stochastic integrals. In sections 4 and 5 we study the forward and backward Taylor expansion, 
respectively. We note that in the stochastic case, the temporal direction of the expansion does 
affect the outcome. Finally in section 6 we try to apply the Taylor expansion to the solution of a 
class of nonlinear SPDEs, and propose a possible definition of the stochastic viscosity solution in 
this case. 

2 Preliminaries and Statement of the Main Theorem 

Throughout this paper we denote (0, T, P) to be a complete probability space on which is defined 
an ^-dimensional Brownian motion B = {Bt)t>o- Let F B = {J~ B }t>o be the natural filtration 
generated by B, augmented by the P-null sets of T\ and let T B = F^. We denote A4 B T to be 
the set of all F fi -stopping times r such that < r < T, P-a.s., where T > is some fixed time 
horizon; and denote -M-q^ to be all F B -stopping times that are almost surely finite. 
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In what follows we write E (also E 1; ■ ■ ■ ) for a generic Euclidean space, whose inner products 
and norms will be denoted as the same ones (•, •) and | • |, respectively; and write B for a generic 
Banach space with norm || • ||. Moreover, we shall denote Q C T B to be a sub-cr-field of T B , and 
for any x G R d and constant r > we denote B r (x) to be the closed ball with center x and radius 
r. Furthermore, the following spaces of functions will be frequently used in the sequel. We denote 

• L P (Q;E) to be all E-valued, ^-measurable random variables £, with < oo. Further, 
we shall denote L^°(£;E) = n p> iLP(£;E). 

• L q (F B , [0, T]; B) to be all B-valued, F B -progressively measurable processes ip, such that 
E Jq \\^t\\ q dt < oo. In particular, q = stands for all B-valued, F B -progressively measurable 
processes; and q = oo denotes all processes in L°(F B , [0,T];B) that are uniformly bounded. 

• C k ' e ([0,T} x E;Ei) to be the space of all Ei-valued functions defined on [0,T] x E which 
are fc-times continuously differentiable in t G [0, T] and £-times continuously differentiable 
in x G E. 

• C fe M ([0,T] x E;Ei), q M ([0,T] x E;Ei), and Cp'\[0,T] x E;Ei), etc. to be the subspace of 
C k ' e ([0, T] x E;Ei), where the subscript "6" means all functions and their partial derivatives 
are uniformly bounded; "/" means all functions are of at most linear growth; and "p" means 
all functions and their partial derivatives are of at most polynomial growth. The subspaces 
with the combined subscripts of b, /, and p are defined in an obvious way. 

• for any sub-CT-field Q C jf, C k > e {G, [0, T] xE; E x ) (resp. C fe M (£, [0, T] xE; Ei), C$' e (G, [0, T] x 
E; Ei)) to be the space of all C k > e ([0,T} x E;Ei) (resp. C^\[0,T] x E;Ei), C^([0,T] x 
E; Ei))-valued random variables that are Q ® B([0, T] x E)-measurable; 

• C^(F B ,[0,T] x E;Ei) (resp. C b M (F B ,[0,T] x E;Ei), Cp' e (F B , [0, T] x E;Ei)) to be the 
space of all random fields 99 G C k ' i (J rB , [0, T] x E;Ei) (resp. C k /{T^ , [0, T] x E;Ei), 
Cp'^F^, [0,T] x E;Ei)), such that for fixed x G E, the mapping (i,w) 1— > (p(t,x,u) is 
F B -progressively measurable, and for P-a.e. oj, the mapping <^(-, -,o;) G C k,e ([0,T] x E;Ei) 
(resp. C fe M ([0,T] x E;Ei), C p M ([0,T] x E;Ei)). 

If Ei = R, we shall drop Ei from the notation (e.g., C fc '^([0,T] x E), and so on); and we write 
C°'°([0,T] x E;Ei) = C([0,T] x E;Ei), and C°'°(F B , [0,T] x E) = C(F B , [0,T] x E),..., etc., to 
simplify notation. Finally, for (t, x, y) G [0, TjxM^xR, we denote D = D x = V X = (^, • • • , ^-) , 
L> 2 = D xx = (d x . Xj )fj =1 , D y = J^, and Z?t = ^. The meaning of D yy , etc., should be clear. 
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Finally, since the random fields that we are interested in are always of the form of (jl.3p . 
which is an SPDE already, the following definition in [5] is useful. Consider the fully nonlinear 
second-order SPDE: 



t rt 

2„ 



u(t, x) = uo(x) + / f(s,x,(u,Du,Du)(s,x))ds + / g(t,x,(u,Du)(s,x))dB s , (2.1) 
Jo Jo 

where (t, x) G [0, T] x R d , and / and g are functions with appropriate dimensions. 

Definition 2.1 A random field u = {u(t,x,uj) : (t,x,uj) G [0,T] x R d x Q} is called a "regular" 
solution to SPDE (OOP if 

(i) ueC°' 2 (F B ;[0,T} xR d ); 

(ii) u is an Ito-type random field with the form 

u(t,x) = u (x) + [ ui(s,x)ds+ [ u 2 {s,x)dB s , (t, x) G [0, T] x R d , 



Jo Jo 

where 

ui(t, x) = f(t, x, (u, Du, D 2 u)(t, x)), U2(t, x) = g(t, x, (u, Du)(t, x)), 
for all (t,x) G [0,T] x R d , P-a.s. U 

In this paper we will consider a special type of "smoothness" for a random field, defined through 
what we shall call the "n-fold differentiability" below. Such a characterization of differentiability 
is mainly motivated by the li stochastic characteristics" for stochastic PDEs (cf. e.g., |12j). which 
often take the form of a system of first order Stochastic PDEs. The idea is actually quite similar 
to the well-known transformation from a higher order ordinary differential equation (ODE) to a 
first order system of ODEs. 

To begin with, we recall that for any multi-index j = (ji, • • • ,jd)i hs "length" is defined by 
lil = Y2t=iJk- We have the following definition. 

Definition 2.2 A random field ( G C°' n (F B , [0,T] x R d ) is called "n-fold" differentiable in the 
spatial variable x if there exist n smooth random fields Q € C°' n (F B , [0, T] xR d ; R di ), 2 < i < n+l, 
with d\ = d and di G N, 2 < % < n + 1, and the functions F i: G{ : R d x -4 R d , i = !,■■■ ,n, 

such that, denoting £\ = the following properties are satisfied: 

(Tl) FuGiZCfi,, i = !,-■■ ,n; 

(T2) Forl<i<n, it holds that 

Ci(t,y) = &M + f F i (yXi+i(s,y))ds+ I Gi{y,Ci + i{s,y))dB s , (2.2) 
Jo Jo 

for all (t,y) G [0, T] x R d , with 0,o G C 2 (R d ;R d <), l<i<n. 
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(T3) For any 1 < i < n + 1, m > 1, and multi-index j = (ji, ■■■,jd) with \j\ = j\ + ... + < n, it 
holds that 

sup{\DiCi(t,y)\,t e [0,T],\y\ < m} e L^(Q,T,P). 

We shall call Ci, i = 2, • • • , n + 1 the "generalized derivatives" of C = Ci? with "coefficients" 
(Fi,Gi), i = !,■■■ ,n. 

For notational convenience we will often write F = F\ and G = G\ when there is no danger 
of confusion. We denote the set of all n-fold differentiable random fields by Cp ([0, T] x M. d ). 

Before we state the main result of this paper, let us note again that the main feature of our 
stochastic Taylor expansion is that it can be expanded around any random point (r, £), and that 
the expansion holds almost surely with an exceptional set that is independent of the choice of 
(r, £). But this amounts to saying that the expansion can be performed around any deterministic 
point (r, £) with any (deterministic) increments, outside a uniform exceptional set. In other words, 
the complicated approximating sequences (ifc,^) m (|1 -2|) can be replaced by simple deterministic 
increments (i + h, x + k), for all (h, k) near zero. We should also note that for the purpose of our 
application, in this paper we consider only the Taylor expansion up to the second order, and for 
that purpose the 3-fold differentiability of the random field would suffice. The precise statement 
of our main result is the following theorem. 

Theorem 2.3 Let £ € C '^ (F B , [0, T\ x M d ) be a random field satisfying the standard assump- 
tions (T1)-(T3) with generalized derivatives Q, i = 2,3,4 and coefficients (Fi,Gi), i = 1,2,3. 
Then, for every a 6 (g, 5) and every m 6 N there exist a subset VL a ,m C £1 with P{£l ai m} = 1, 
such that, for all (t,y,uj) £ [0, T] x B m (0) x 0, atTn , we have the following Taylor expansion 

C(i + h, y + k) - C(t, y) = ah + b(B t+h - B t ) + ^(B t+h - B t ) 2 + <p, k) + ~{Xk, k) (2.3) 

+(q, k)(B t+h - B t ) + (\h\ + \k\ 2 ) 3a R atm {t, t + h,y,y + k), 

for all (t + h,y + k) G [0, T] x B m (0). Here, with (F, G) = (F\,G\), one has 

a = F(y,( 2 (t,y))-^(D z G)(y,( 2 (t,y))G 2 (yX3(t,y)), 

b = G(y,( 2 (t,y)), c=(D z G)(y,( 2 (t,y))G 2 (y,(3(t,y)), 

p = D y ((t,y), X = D 2 y ((t,y), (2.4) 

q = (D y G)(y,(2(t,y)) + (D z G)(y,(2(t,y))D y ( 2 (t,y). 

Furthermore, the remainder of Taylor expansion R at m '■ [0,T] 2 x (IR^) 2 x Q 1— > M. is a measurable 
random field such that 

R a>m = sup _ \R a)m (t,s,y,z)\ eL™(n,F,P). (2.5) 

t,se[Q,T\; y,z&B m {0) 
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Proof. Since the proof of this theorem is quite lengthy and technical, we shall split it into 
several cases and carry it out in the following sections. ■ 



Remark 2.4 (i) It is worth noting that the "universal" estimate for the remainder R a , m is the 
main reason why the Taylor expansion can now hold around any random space-time point. In 
other words, the points (t,y) and (h,k) in Theorem 12.31 can be replaced by any (r, £), (cr,rj) € 
L°(F B ; [0,T]) x L°(F B ;R d ), such that (r + 0-^ + 77) € [0,7] xB m (0), P-a.s. on Q, except in that 
case the reminder should read 

R a ,m = Ra,m(r,T + cr, £ + rj) , where \R a , m \ <R a , m £ L°°(F B ;P). 

In what follows we denote R a , m to be a generic term satisfying (j2.5j) . which is allowed to vary 
from line to line. 

(ii) Prom the expressions (12, 4ft it is clear that the drift term F(= F±) appears only in the 
coefficient a, and it is in its original form. In fact, as we shall see in the proof, the exact form 
of F(t,y)(= F(y,(2(t,y))) is not important at all. We could simply change it to F(t,y) and the 
results remains the same. 

(iii) Theorem 12.31 remains true if the assumption (Tl) is replaced by the weaker assumption: 
(Tl') Ii, (I/: Q„ ? : L<J2 3. 

However, we shall not pursue this generality in this paper due to the length of the paper. ■ 

To end this section we give an example, which more or less motivated our study. 

Example 2.5 Let u = {u(t,x)} be a regular solution of the SPDE (12. ip . and for simplicity we 
assume that both / and g are "time- homogeneous" (i.e., they are independent of the variable t), 
and g is independent of x as well. We define £ = £1 = u, and Ci+i = (C«) F>Q, D 2 Q), 1 < i < 3. 
Then, assuming that the coefficients / and g are sufficiently smooth and their derivatives of all 
order are bounded, one can show that (T1)-(T3) are satisfied. Furthermore, applying Theorem 
12.31 we see that on some subset of £1 with P(£l) = 1, independent of the expansion point 



S 



(t,x) G [0,T] x R d , the stochastic Taylor expansion ()2.3p holds for u, with 

a = /(x,(u,Du,Z)u 2 )(t,x))-i{( 5j D u ( ? )((n, J Du)(t,x)) 

+A.s((u, £>u)(i, x))(D p g((u, Du){t, x)),Du{t, x)) 

+ (D 2 u(t, x)D p g{{u, Du)(t, x)),D p g({u, Du)(t, x))}} 
b = g((u,Du)(t,x)) 

c = (gD u g)((u,Du)(t,x))+D u g{(u,Du)(t,x)){D p g{(u,Du)(t,x)),Du(t,x)) (2.6) 

+(D 2 u(t, x)D p g((u, Du)(t, x),D p g((u, Du)(t, x))) 
p = Du(t,x); X = D 2 u(t,x) 

q = D u g((u,Du)(t,x)Du(t,x) + D 2 u(t,x)D p g((u,Du)(t,x)), 

for all h G R and k G M. d such that t + /i G [0,T] and |&;| < m, respectively. Here, the coefficients 
F2,G2,F^, and G3 are obtained by differentiating f)2. 1 j) with respect to x. 

The explicit expressions of the pathwise Taylor expansion in (|2.6p will be important for our 
study of the stochastic characteristics for SPDE (|2. 1 j) . which will be discussed in the last section 
of this paper. ■ 



3 Some fundamental estimates 

Before we prove the theorem, let us first introduce the so-called Wick-square of the Brownian 
motion, which is originated in the Wiener homogeneous chaos expansion (cf., e.g., [9], |19j). For 
any < s < t, we define the Wick-square of B t — B s to be 

(B t -B s f 2 := (B t - B s ) 2 - \t — s\. 

Moreover, let L^ B (Q; W 2 °'/ oc ([0, T] x R d )) denote the space of all random fields / G L°(F B , [0, T] x 
R rf ) such that /(w, •, •) is P-a.s. an element of 

W^ oc ([0,T] x R d ) = {uG L 2 oc ([0,T] x R d ) : D x u G L 2 oc ([0,T] x R d )}, 

where in this case D x u denotes the weak partial derivative of u with respect to x. 

We begin by proving an important estimate for multiple Stochastic integrals. We consider the 
multiple stochastic integral defined recursively as follows. Let fj G L^ B (Q; / oc ([0j T] x R d )), 
for 1 < j < N, N G N. Let X^ s (x) = 1, and for j = 1, ■ ■ ■ , N, 

Xi s (x) = j° f 3 (s 3 ,x)xi^(x)dB Sj , se[t,T}. (3.1) 
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It is also useful to define the similar multiple integrals: for k > 0, let X t £ 1 (x) = 1, and for 
I > k, 



x k A{x) 



j fi{u h x)X^l l 1 (x)dB Ul , se [t,T]. 



(3.2) 



Then it is clear that X t g = Xf s , for j > 1. We have the following regularity result. 



Proposition 3.1 Let JVgN be given. Assume that fj G L° s (O; W 2 ' ioc ([0, T] x 
satisfy that 



,l<j<N, 



Cm,p(/j) — E 



x\<m J 



T 



\Difj(t,x)\ p dtdx 



< oo, 



(3.3) 



for all m G N, p > 1, and i = . . . , id) with \i\ = i\ + . . . + < 1. Then it holds that 



C/3, m = sup<^ , _U , 0<t <s<T, \x\<m \ €L°°(Q,T,P) 



(3.4) 



/or all (3 G (0, f ), m G N. 



Proof. The proof is based on the Kolmogorov continuity criterion for random fields (cf. e.g., 
Theorem 1.2.1]), combined with an induction argument. We note that throughout the proof 
we shall use the notations C m , C mjP , CM,k, etc., to represent the generic constants depending only 
on /j's, T, and the parameters in their subscripts, which are allowed to vary from line to line. 

We begin with the case N = 1. In this case we write fi = f and X\ = X for simplicity. Let 
m G N, 7 > 2, and < t < s < T. Applying Sobolev's embedding theorem, Holder Inequality, 
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and noting that X^ s = 1, we have 



< 



< 



< 



< 



E 



sup |A" M (x)[ 7 

\x\<m 



< C m E 



£ / \DiX t)S {x)\ d+1 dx 
^ sl J\x\<m 



d+l 



|J|<1 



si <m 



Dlf(r,x)dB r 



7(d+l)' 



l 

d+l 



h'[<l Wr^ 71 ^ Jt 



\D J x f(r, x)\ 2 dr 



7(d+l) 
2 



1 

d+l 



bi<i 



7(d+l) 



(s-i) 7 ^!) 2 ^ |L)i/(r,x)| 7(d+1) dr 



7 (d+l) 



c^E (7 



{s-t) 1 ^- 1 I \Dif(r,x)\^ d+1 Ur 



dx 



< C m (s-t)--^ £ I -E 
b1<i V 



|i|<m ^0 



i 

d+l 



Consequently, 



sup |X tiS (x)|' 



cc <m 



<C m , 7 (s-t)2 <J+i. 



Endowing the space C(B m ) of continuous functions defined on the closed m-ball B rn := {x E M. d : 
\x\ < m} with the sup-norm, and considering {x h-> Xo, s (a;); |sc| < w-} s e[o,T] as a C(i? m )-valued 
process, we see that 



sup \X t>s (x)\ = \X ,s(-) - X ,t(- 

\x\<m 



\C(Bm)- 



Applying the Kolmogorov continuity criterion, we conclude that C/3,m € (f2, J 7 , P) for all f3 £ 
0, | — 7 |^ 2 1 ) 1 • Since we can choose 7 > 2 arbitrarily large, f|3.4[) follows. 

We now prove the inductional step. That is, we assume that (|3.4p is true for N — 1, and show 
that it is also true for N. To do this, we shall adapt the proof of the Kolmogorov continuity 
criterion given in [23] to our framework. Let < t < s < T, 7>2, and m € N. First note that 
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by a simple application of Burkholder-Davis-Gundy and Holder inequalities we have 



E{ sup \X*{x)\A = E{ sup f f N (s N ,x)X^ SN 1 (x)dB SN \ P \ 
1 re[t,s] } 1 re[t,s] Jt 1 } 

< C p E{[J°\f N (s N ,x)\ 2 \Xl-\x)\ 2 d SN \ Pl2 } 



(3.5) 



< C p E{\\f N (.,x)\\l 2M sup \X t N s -\xW] 



< a 



' P {E\\fN(; 



*)IISm)} 1/A { e { s T J x ^)\^ 



By simply iterating the above argument we obtain that 



E{ sup \X»{x)T} 

L rE[t,s] ' 



N 



:,pN 



pN i AT-2 



< ^H/atO^)^ [t , ■ EU^x)]^ M -E{ sup Xf- 1 ( 2 ;)|^} J (3.(5) 



< 



< 



N 



k=l 



Furthermore, with a similar argument as above we can also show that for all M, k € N, p > 1, 
and multi-index i satisfying \i\ = 1, 



e{ sup |M t ^(a;)| fe 



< El 



D iM f M ( / • • • ( / D' J\>IIL) ■ ■ ■ dB SM1 J dB SM 



< c M>k II 7 /r •/>> 

|ii|,...,|i*f[<li=l 

M 

^ n^[p ij /,-fe^))iifA[*, 

|i 1 |,...,|i M |<n=i 



A/fc(p-l) 
< C M>fe (s-t) 2p 



1 



Here, we applied the Holder inequality in the last step above. Consequently, if p < then the 



assumption (|3.3p implies that, for all m > 0, 



L 



E 



|x|<m 



sup ^X^ixt 

t<r<s 



Affc(p-l) 

dx < C mtM>k (s - t) 2 p 



(3.7) 
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We can then conclude that, for all < t < s < T, 7 > 2, 



E 



snp 

\x\<m 



< C m E 



< C m E 



£ / \DiX t N s (x)\ d+1 dx 

mi J\x\<m 



d+l 



E 

vN.lil<i 



x <m 



DiX»- l (x)Bif N (r,x)dB r 



y(d+l) 



1 

d+l 



< G 



< 



|x|<m 



v|i|,li[<l" 



xl <m 



sup 

t<r<s 



dr 



1 

d+l 



Consequently, by Holder inequality and (]3.7[) with p = -s- = (iV — l)"f(d + 1), we have 



< C n 



sup ixf s (x)r 

|a;|<m 



N,lil<i 



x|<m L i<r<s 



sup 



x|/ E 

\x\<m 



< 



Gn, 7 £ ((*-*)* 
|i|<l,i=l,2 V 



2 (jV-l) 7 (d+l)-§ x 



dx 



d+l 



x I / E 

\x\ <m 



7(J+D-i . _ \7(<i+l) 

(*-*) l{d+l) \\ D lfj(-^)Wl^)([t,s]) 



2 \ d+l 

dx 1 



d+l 



E * 



W<1, j=l,2 



|z|<m 



in' f.f. T ^|i 2 7(rf+i) j 

l^a-M 5 -^y 1 II ^2T(d+l) ([t >s ]) 



In other words we obtained that 



E 



sup 

\x\<m 



< Gn, 7 (s - t)^-^ , 0<t<S<T. 
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Next, recall the multiple integrals X^' l s {x) denned by ()3.2p . For < r\ < r 2 < r 3 < T, and 
1 < n o < N, define 

Y ri,r2,r3,no 0*0' X T2 ,r3 0*0-^ri,r2 0*0 > (3-9) 

An easy calculation shows that 



„ r3 
/ n 



r2 „ r3 

x l[ N UN 1 {x)f N {u N ,x)dB UN + / X^~ 1 (x)f N (u N ,x)dB UN 

n J r 2 

JV-2 

= A^ 2 (x) + ^ ^ r2 , r3 , no (^) + X^ix) (3.10) 

n =0 

To simplify the further discussion, we now assume without loss of generality that T = 1. Let 
D n be the set of all tf := 2~ n i, for some < i < 2 n . The set D = [j n>l D n is then the set of all 
dyadic numbers in [0, 1]. Now, let n G N and 0<i<s<lbe arbitrary dyadic numbers such 
that s — t < 2~ n . Our aim is to estimate Xt tS (x) uniformly in t,s G D. To this end, we notice 
that there exists some k > n such that t, s belong to D^. Moreover, denoting 

Sj = sup{r G Dj, r < s}, tj = sup{r G Dj, r < t}, n < j < k, 
we have s n < s n+ i < . . . < s^ = s, t n < t n+ i < . . . < = t. Therefore (see also Figure [TJ 

fc-l ( N-2 } 

*i» = 4:1 0*0 + E + £ Y ts 3 ^nM , (3.H) 

j=n (. n =0 ) 

fc-l ( N-2 *\ 

X$(x) = Xlf(x) + Y, X];l +1 0*0 + E Y Ct 3+ us,n (x) • (3-12) 

j=n I n =0 ) 

In order to estimate Y£ j8 . j8 . +uno (x) and Y t f jt . +1)S}nQ (x), we notice that 

s j+1 - s 3 G {0, 2-^' +1 )}, t j+1 - tj G {0, 2-^' +1 )}, n<j<k, 

and thus 



(3.13) 



f < Sj - t n = (Sj - s n ) + (s n -t n )<2- 2~ n ; 
\ < a - t j+ i = (s-t) + (t- t j+1 ) < 2 • 2~ n , 

for n < j < k. Moreover, recall that for every a G (0, ^), there exists a 9 a G L^°(f2, J 7 , P) such 
that for all /, /' G N with </'-/< 2V - 2, 

P { X l ^ +h < 6 a h {l '- l+1)a , < i <t + h < T, |x| < m} = 1, 
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Y t m ,s l5 s 1+1 (y) 



s l s l+l s ~ s k 



Figure 1: X Sjt = X]f and lt, SljS2 = ^ Sl)S2) n for iV = 2 and n = 1 



thanks to the inductional hypothesis. Thus, we have by f)3. 13|) 



N-2 

n , Sj ,s j+1 ,n { x )\ 

no=0 



N-2 



< <*\s j+ i - S ,| ( " 0+1)Q • e a \ Sj - t n \( N - n °-^ a 

no=0 
JV-2 

< ^2 9 2 2~ ( - no+1 ^ j+1 ^ a 2~ ( - N ~ no ~ 1 ^ n ~ 1 ' )a 

no=0 

N-2 

< 2 N ^ 02 2~((n +l)j+(N-n -l)n)a 

no=0 



(3.14) 



Similarly, we obtain 



N-2 

E 

no=0 



N-2 



\Y N (x)\ <2 N fl 2 2-(( n o+ 1 )j'+( iv - n o-l)n)a 



/ j "a 
no=0 



(3.15) 



Next, we define Kj = sup{|X j j (x)\, \x\ < m, < i < 2 J — 1}, n < j < k. Then, combining 
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(13TT|) . ([3~T2]) . (|3~T5t and (|3~Tlt we deduce 



\X t N s (x)\ 



k-l 



N-2 



N-2 



< Sn w + E x 5U+» - x Zt j+1 w + E y * 



N 

n ,Sj,s j+1 ,n 



(x) E ^i i+ i,s,n ( a; ) 



j=n 



no=0 



n =0 



fc-1 



JV-2 



< k„ + e 2^ + 2^E^ 2 " 



((n +l)i+(AT-no-l)n)a 



no=0 



N-2 



< ^ 3^. + 2^+! ^ 6> 2 i 2~ ((no+1)i+(Ar ~"' ~ 1)?l)a 

n =0 / 



j=n 



Now let (iV - l)a < /3 < Na < f . Define 

C^^supj 1 ^^ 1 , s,t e D, 2-(™ +1 )< S -t<2-", |x|<m|, 



A 



and Mp im = sup n > 1 C« m - Then, 



Mp >m < sup (2 (n+1)/3 sup{|X 4 iV s (x)|, s,teD, 0<s-t<2~ n , \x\ < 



n>l 



in 



N-2 



< sup 2(" +i )^ e 3K i + 2N+1 E 



((no+l)j+(Ar-no-l)n)a 



n>l 



j=n 
oo / 



no=0 



JV-2 



< sup2 /3+Ar+2 E 2 nj8 ifj + E ^2" (no+1)Qj+(/3 " (7V ~ n °~ 1)a)Tl 

n — 1 j=n \ rao=0 / 

oo / JV-2 > 

< sup2 /3+Ar+2 E ( 2 nl3 Kj + E 2 2 _(no+1)QJ ' +( ^~ (w ~ 1)a)n+rm ° Q 

j=n \ no=0 / 



n>l 



oo 

< 2? +N+2 E (^ l3 K J + (AT - 1)^2^-^') 
i=o 

= 2^ N+2 (y^K j+ {N ~^° l 



Noting that by ([M 



2^-1 



< E E 



i=0 



sup |Xi j (x)| 7 



|x|<m 



<2^ m , 7 2^f^)=C m , 7 2-^f^ 2 ), 



(3.16) 
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we thus obtain for 7 > 



N-2/3 



that 



E 



< 2 



P+N+2 



+ 



N -I 



J=0 



2 Is 



1 _ 2-Wol-P) 

N - 1 



+ 



m,7 



+ 



1 _ 2 -(Na-P) 

N-l 



E[9%]) 

{e[oV]Y 



< 00. 



Consequently, Mp >m G L^°(f2, J 7 , P), for all /3 G (0,Na). But since a € (0, |) is arbitrary, we can 
extend the result to /3 G (0, N/2). Finally, note that the definitions of Mg iT71 and {Cfl m }i as wen 
as the continuity of the mapping (t, s) i-> X( )S (x), imply that 

'\X t N s (x)\ 



C/3,m = SUp 



0<s-t<2~\ \x\ <m \ £L™(fl,7,P), 



, l^-t| /3 

for all /3 G (0, iV/2) and m G N. The proposition then follows from the recursive relation lOD .H 

The following corollary can be easily obtained by adapting the proof of Proposition 13.11 in an 
obvious manner. 

Corollary 3.2 The statement of Proposition [Q1 remains valid if, for 1 < i < N, dB Si is replaced 
by dA\. in \3. where A 1 is either the Brownian motion B or A\ = s, s G [0, T\. Moreover, 
13. 4\ ) holds whenever (3 G (0,l\ + , m £ ¥S , where l\ is the number of i, 1 < i < N , for which 
A\ = s, s G [0,T], and l 2 = N - h. 



4 Forward Taylor expansion 

In this and next section we shall provide a complete proof of Theorem 12.31 The results in each 
section, however, can also be applied independently, and therefore are of interest in their own 
right. We should note that, unlike the usual Taylor expansion, in the stochastic case the direc- 
tion of the time increment makes significant difference in the argument, due to the "progressive 
measurability" of the random fields. We will thus separate the two cases: in this section we study 
the forward expansion, and leave the backward case to the next section. 

4.1 Forward temporal expansion 

We begin with the forward "temporal" expansion, that is, only the time variable has the increment. 
Let us first introduce the following extra notations. For each T > we define the 2-dimensional 
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simplex 

A [0)T] = {(t, s) G [0, T] x R+ : < t < s < T}. (4.1) 

For any real- valued measurable functional 6 denned on x Amy] x IR d+1 and m G N, we denote 

9 m = sup{\9(t,s,x,z)\, (t,s) G A [0jT ], x G R d , z G R with |z| < m}. (4.2) 

Furthermore, by a slight abuse of notations, in what follows we shall denote R a to be any mea- 
surable functional 9 a , indexed by a € (|, ^), such that (9 a ) m G L°^{Vt,F,P) for all m G N, and 
again, it may vary from line to line. 

Our main result of this section is the following Stochastic forward temporal Taylor expansion. 



Proposition 4.1 Assume that ( G C°' (3) (F S , [0, T] x R d ) satisfies (Tl)-(TS). Then, for all 
G IN and a G Q, |) t/iere exists a subset fl C ft such that P(£l) = 1, and that on f2, for all 



m 



0<t<t + h<T,x£ Mr, the following expansion holds 

C(t + h,x)- C(t, x) = ah + b{B t+h - B t ) + |(P m - P t ) 2 + h 1+a R a (t, t + h,x), (4.3) 

a = F{x,( 2 (t,x))-^(D z G)(x,( 2 (t lX ))G 2 (x,t 3 (t,x)), 

b = G(x,C 2 (t,x)), (4.4) 

c = (^G)(x,C 2 (t,x))G 2 (x,C3(t,x)). 

Proof. Let /i > be such that 0<t<t + /i<T. For any x G R d , we write 

rt+h pt+h . 

C(t + h,x) -((t,x) = F(x,(2(s,x))ds+ G(x,( 2 (s,x))dB s = I 1 + 1 2 , (4.5) 

where P = P(t, h, x), i = 1,2, are the two integrals. We shall study their expansions separately. 

We begin by I 1 . The argument is very similar to that of [3], we provide a sketch for complete- 
ness. Let 

H\x, z 2 , z 3 ) = (D z F)(x, z 2 )F 2 (x, z 3 ) + hi[G 2 (G 2 )*(x, z 3 )D 2 z F(x, z 2 )} 
for (x, z 2 , z 3 ) G R d x R^ 2 x R ds . Then, applying Ito's formula, we have 
F{x,(2{s,x))-F{x,( 2 (t,x)) = [ P 1 (x,C 2 (r,x),C 3 (r,x))dr 



+ / ((D z F)(x,( 2 (r,x)),G 2 (x,( 3 (r,x)))dB r , (4.6) 
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for < t < s < T and x G R d , P-a.s. We now show that there is a universal subset lycfl with 
P(yi') = 1, on which (|4.6|) holds for all < t < s < T. But for this it suffices to prove the random 
fields 

I" H 1 (x,C 2 (r,x),C 3 (r,x))dr, [ {(D z F)(x, C 2 (r, x)), G 2 {x, ( 3 (r, x)) ) dB r 
Jt Jt 

have continuous versions. To see this, we try to make use of the Kolmogorov continuity criterion. 

For notational simplicity let us denote 

T(t,x) = (D z F){x,C2(t,x))G 2 (x,&(t,x)), (t,x) G [0, T] x R d . 

Then, we can deduce from (T3) and the Burkholder-Davis-Gundy inequality that, for each k > 2, 



E 



T(r,x)dB r - / T(r,x')dB. 



< CUE 



|(r(r, x)\ 2 dr 



+ E 



\T(r,x)\ 2 dr 2 + Ei I \T(r, x) - T(r, x')\ 2 dr 



< C k (\t' + \ s ' -a\%) El sup \T(r,x)\ k \ + C k \x-x'\ k E\ sup \D x T(r,x)\ k \ 
\ ' \ re[o,T], J [ re[o,T], J 



< Cfc ( \t' - t\2 + Is' - S\2 + [x - X 



Ilk 



for < s < t < T, < s' < t' < T, \x\, \x'\ < m, and k G N. Hence, the Kolmogorov continuity 
criterion renders that the random field | T(r,x)dB r , < t < s < T, |x[ < m| possesses a 
version that is continuous in (t, s, x). A similar estimate allows to prove that also the first integral 
in (|4.6p admits a version continuous in (t,s,x) as well. Hence, we conclude that on some f2' C fi 
of full probability measure the relation (|4.6p holds for all < t < s < T and < m. 

Consequently, writing the integral I 1 as 

/t+ft rs 
^ H\x,C 2 (r,x),C 3 (r,x))drds 

+ [ + f (D z F)(x,(2(r,x))G 2 (x,(3(r,x))dB r ds, 
Jt Jt 

and using the conclusion above we see that as a function of (t,h,x), I 1 is jointly continuous for 
allO<i<i + /i<T and x G M d , over the universal set fi'. It then follows from Corollary 13.21 
that there exists an 9," C Q' with P(0") = 1 such that for all < t < t + h < T and x G M ', it 
holds that 

I 1 ^, h, x) = F(x, ( 2 (t, x)h + h 1+a R a (t, t + h, x). (4.7) 
In what follows we will not distinguish Vt" from 0,'. 
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We now turn our attention to I 2 . Again, we begin by denoting, for (x, z 2 , £3) £ R rf x R^ 2 x M^ 3 , 
H 2 (x,z 2 ,z 3 ) = (D z G)(x,z 2 )F 2 (x,z 3 ) + ^tr[G 2 (G 2 y(x,z 3 )D 2 G(x,z 2 )}. (4.8) 
Then, for every x € R d , we again apply Ito's formula to get, for < t < s < T, and P-a.s., 

G(x,(2(s,x))-G(x,(2(t,x)) = J* H 2 (x,Ur,x),Ur : x))dr (4.9) 

+ [ S {D z G){x,C, 2 (r,x))G 2 {x,Ur,x))dB r . 



Using the similar arguments as before we can find another universal subset, still denoted by 
n'Cfl with P(O') = 1, on which flUD holds for all < t < s < T. 
Next, using (|4.9p it is easy to see that I 2 can be written as 

I 2 (t,h,x) = G(x,£2(t,x))(B t+h -B t )+ / H 2 (x,(2(r,x),(3(r,x))drdB s 

Jt Jt 

+ / / (D z G)(x,(2(r,x))G 2 (x,(3(r,x))dB r dB s 



t Jt 



= G(x,( 2 (t,x))(B t+h -B t ) + I 2 ' 1 {t,h,x) + I 2 ' 2 (t,h,x), (4.10) 

and we can claim as before that (|4.10p holds for all < t < t + h < T, \x\ < m, over a universal 
subset of full probability measure, again denoted by Q' . 

We now analyze I 2 ' 1 and I 2 ' 2 separately. Using integration by parts we see that 

I 2 ' 1 = (B t+h - B t )£ +h H 2 (x,(2(r,x),( 3 (r,x))dr \B r -B t )H(x, ( 2 (r, x), ( 3 (r, x))dr. 

Then, following the argument developed in the previous part, we can show that the equality holds 
for all < t < t + h < T, x G R d , over anfi'Cfi with P(Q') = 1. Note that 

sup {h- a \B t+h -B t \}eL°°(n,T,P), (4.11) 

0<t<t+h<T 

it follows that over £1' , 

I 2 ' 1 ^, h, x) = h 1+a R a (t, t + h,x). (4.12) 
The estimate for J 2 ' 2 is slightly more involved. For notational simplicity let us define 



K(x,z 2 ,z 3 ) = (D z G)(x,z 2 )G 2 (x,z 3 ), (x, z 2 , z 3 ) e R x R 2 x 



F(x,z 3 ,z 4 ) = (F 2 (x,z 3 ),F 3 (x,z 4 )), (x, z 3 , Z4) e R x R 3 x 



A 



G(x,z 3 ,z 4 ) = (G 2 (x,z 3 ),G 3 (x,z 4 )), 

H(x,z 2 ,z 3 ) = D ?2 K(x,z 2 )F(x,%) + ~tr[GG%x,%)DlK{x,z 2 )], 
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where z 2 = (#2,23) an d % = (2:3,214). Moreover, we denote Q(s,x) = (Q(s,x),Q + x(s,x)), i = 2,3. 
Then, applying Ito's formula we have, 



K(x,(2(s,x)) ~ K(x, C2(t,x)) = J H(x,C2(r,x),(3{r,x))dr 

+ J" Dz 2 K(x, C 2 (r, x))G(x, C 3 (r, x))dB r . 
Again, we assume that the equality holds for all < s < t < T and x G M. d , on Therefore 

'■t+h rs 



rt+n rs 

I 2 ' 2 (t,h,x) = / K(x,( 2 (r,x))dB r dB s 

Jt Jt 



H Jt 

^(D z G)(x,( 2 (t,x))G 2 (x,( 3 (t,x))(B t+h - B t f 2 (4.13) 

rt+h rs rv 



It Jt Jt 

Moreover, Proposition 13.11 implies that 

'■t+h rs 



rt+n rs rr 

+ / / / H(x,( 2 (v,x),(3(v,x))dvdB r dB s 
Jt Jt Jt 

/•t+h rs rv 

+ / / / D^K(x,( 2 (v,x))G(x,( 3 (v,x))dB v dB r dB s 
Jt Jt Jt 



rt+n rs rr 

k-(i+a) / / / D % K(x,( 2 (v,x))G(x,C3(v,x))dB v dB r dB s = R a (t,t + h,x); 
Jt Jt Jt 



't Jt Jt 
and Corollary 13.21 implies that 

rt+h rs rr 

h -(i+2a) / / H(x,( 2 (v,x)X 3 (v,x))dvdB r dB s = R a (t,t + h,x). 

Jt Jt Jt 

Consequently, we see that, over £1', (|4.13p becomes 

I 2 > 2 (t, h, x) = ^{D z G){x, ( 2 (t, x))G 2 (x, C 3 (t, x))(B t+h - B t f 2 + h 1+a R a (t, t + h,x). (4.14) 

Finally, plugging (@2J), KTU\i . <^A2\f . and (jlTUD into g3D, we obtain and flH]), with a 
universal exceptional null set, proving the proposition. ■ 



4.2 Forward temporal-spatial Taylor expansion 

Based on the forward temporal Taylor expansion Proposition 14.11 we now add the spatial incre- 
ment. Our main result of this section is the following proposition. 

Proposition 4.2 Assume that ( G C°' (3) (F B , [0,T],M d ) satisfying (Tl)-(TS). Then, for all m G 
IN and a G (|, ^) there exists some subset f2 c such that P(&) = 1, and that on (}, for all 
< t < t + h < T and x, k G R d , 

C(t + h,x + k)- C(t, x) = ah + b(B t+h - B t ) + - B t ) 2 + (p,k) +^(Xk, k) 

+ ( q, k)(B t+h - B t ) + (h+ \k\ 2 ) 3a R a (t, t + h,x,x + k), (4.15) 
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where 

a = F(x,t2(t,x))-^{D z G){x 1 (2{t,x))G 2 {x,(3{t,x)); 

b = G(x,(2(t,x)); c=(D z G)(x,C, 2 {t,x))G 2 {x,Qz{t,x))- (4.16) 

p = (D x ()(t,x); X = D 2 x C(t,x); 

q = (D x G)(x,( 2 (t,x)) + (D z G)(x,( 2 (t,x))D x ( 2 (t,x). 

Proof. First let us write 

C(t + h,x + k)- C(t, x) = [C(t + h, x + k) - C(t + h, x)] + [((t + h,x)- CO, ar)], 

where the second [• • • ] above is the forward temporal expansion studied in the previous subsection. 
In light of Proposition 14.11 we need only prove the following stochastic spatial Taylor expansion: 

C(t + h,x + k)-((t + h, x) = D x ((t, x)k + -(Dl((t, x)k, k) + {(D x G)(x, ( 2 (t, x)) 

+ (D z G)(x,( 2 (t,x))D x ( 2 (t,x)}k(B t+h - B t ) (4.17) 
+ (h+\k\ 2 f a R 1 a (t,t + h,x,x + k). 

To this end, we first apply the standard Taylor expansion and use the assumption (T3) to get 

C(i + h, x + k) - C(t + h, x) = D x ((t + h,x)k + -(Dl((t + h,x)k,k) 

+\k\ 2, R 1 a (t,t + h,x,x + k), (4.18) 

for all < t < t + h < T, x, k E R d , P-a.s. Next, differentiating the equation for £ = Ci in ([23]) 
we have 

D x ((t,x) = D x ( (x)+ [ {(D x F)(x,(2(s,x)) + (D z F)(x,(2(s,x))D x ( 2 (s,x)}ds 

Jo 

+ / {(D x G)(x,(2(s,x)) + (D z G)(x,( 2 (s,x))D x ( 2 (s,x))}dB s , 
Jo 

for all t G [0,T], x € M. d , P-a.s. Now, applying Proposition 14.11 (to D x () one can check that 
D x ((t + h,x) 

= D x ((t, x) + {(D x F)(x, ( 2 (t, x)) + (D z F)(x, ( 2 (t, x))D x Q 2 {t, x)}h 

-^{D x [(D z G)(xX2(t,x))]G 2 (x,C3) + D z G(x,C 2 (t,x))D x [G 2 (x,C3)]}h 
+{(D x G)(x,(2(t,x)) + (D z G)(x,( 2 (t,x))D x ( 2 (t,x)}(B t+h - B t ) 
+l{D x [(D z G)(x, ( 2 (t, x))]G 2 (x, C 3 ) + D z G(x, £ 2 (t, x))D x [G 2 (x, ( 3 )}}(B t+h - B t f 2 
+h 1+a R a (t,t + h,x), 
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for all < t < t + h < T, x € R d , P-a.s. Consequently, it follows from ([jXTjl that 

D x ((t + h,x) = D x C(t,x) + {(D s G)(x,C 2 (t,x)) + (D z G)(x,£2(t,x))D x ( 2 (t,x)}(B t+h - B t ) 
+h 2a R a (t,t + h,x), (4.19) 

for all < t < t + h < T, x € M. d , P-a.s. Similarly, one shows that 

Dl((t + h,x) = Dl((t,x) + h a R a (t,t + h,x), < t < t + h < T, x € R d . (4.20) 

Combining KWi . (gUJ), and gSOD , we obtain that 

C(t + h,x + k) -C(t + h,x) 
= D x ((t + h,x)k + ^(Dl((t + h,x)k,k) + |A;| 3 ^(t,t + h,x,x + k) 

= D x ((t,x)k + ~{Dl((t,x)k,k) 

+{(D x G)(x, C 2 (t, x)) + (A,G)(z, C 2 (t, x))^C 2 (t, s)}fc(S t+ h - B t ) 

+ |fc| 3 i?^(t, i + /i, x, x + fc) + h 2a R a (t, t + /i, x)A; + h a (R a (t, t + /i, x)k, k), 

for < t < t + /i < T, x, G R d , P-a.s. Finally, noting that 

[Af + h 2a \k\ + h a \k\ 2 < C m (h 3a + |A:| 3 ) < C m (h + \k\ 2 f a 

for all h £ [0,T], \k\ < m, we see that (14.17P holds, hence the proposition follows. ■ 

5 Backward Taylor expansion 

In this section we treat the backward Taylor expansion, that is, when the temporal increments 
are negative. As a general belief such an expansion would be more difficult than the forward one, 
due to the obvious "adaptedness" issue. But we shall see, with our "pathwise" approach, such 
difficulty is eliminated. We nevertheless would like to separate its proof from the forward case 
because of the slight difference in the arguments. We again take two steps: first the backward 
temporal expansion, and then the mixed time-space expansion. 

5.1 Backward temporal expansion 

We have the following analogy of Proposition 14.11 



23 



Proposition 5.1 Assume that ( G C°'( 3 )(F B , [0,T],M d ) satisfies (Tl)-(TS). Then, for allm£N 
and a € (|) 5) > there exists some subset Cl C £1 with P(Cl) = 1, such that on Q, for allO < t — h < 
t < T , x € M. d , the following expansion holds 



C(t -h,x)- C(t, x) = ah + b{B t _ h - B t ) + -(Bt-. h - B t f + h 1+a R a {t - h, t, x). (5.1) 



where 



a = - lF(x,( 2 (t,x)) -^(D z G)(x,( 2 (t,x))G 2 (x,(3(t,x))\; 

b = G(x,( 2 (t,x)); (5.2) 

c = {D z G)(x,(2(t,x))G 2 (x,£3{t,x). 

Proof. Applying the forward Taylor expansion Proposition 14. 1\ we have 



((t,x)-({t-h,x) 



^F(x,(2(t - h,x)) - ^(D z G)(x,(2(t - h,x))G 2 (x,&(t - M))} h 



+G(xX2(t-h,x))(B t -B t _ h ) 

+^(D z G)(x,( 2 (t - h,x))G 2 (x,( 3 (t - h,x))(B t - B t _ h ) 2 
+h 1+a R a (t- h,t,x), (5.3) 

for all < t - h < t < T, x £ R d , on a full probability set W. Our main task is to replace the 
temporal variable t — h by t. To do this, we first apply Ito's formula to G(x,C 2 (t,x)) to obtain 



G(x,( 2 (t,x)) = G(x,(2{0,x))+ / H(x,( 2 (s,x))ds + L{x,C, 2 {s,x))dB 8 , (5.4) 

Jo Jo 

where ( 2 (t,x) = (( 2 (t,x),( 3 (t,x)), and 

H(x,(z 2 ,z 3 )) = (D z G)(x,z 2 )F 2 (x,z 3 ) + ^t[G 2 (G 2 )*(x,z 3 )(D 2 z G)(x,z 2 )}, 

A (,5.0J 
L(x,(z 2 ,z 3 )) = (D z G)(x,z 2 )G 2 (x,z 3 ). 

Next, we denote z = (22,^3), and define ( 3 (t,x) = (( 3 (t,x),(4,(t,x)), G(t,x) = G{x,(, 2 (t,x)) and 

G 3 (x,(z 3 ,z A )) = (G 2 (x,z 3 ),G 3 (x,z 4 )). 
Applying Proposition 14. 1 1 to (?(•, ■) and using (I5.4p . we deduce that 

G(t,x) -G(t-h,x) = !.H(x,( 2 (t-h,x)) - ^D z L(x,( 2 {t - h,x))G 3 {x,( 3 (t - /i,x))j h 

+L(x,( 2 (t-h,x))(B t -B^ h ) 

+^D z L(x, ( 2 (t - h,x))G 3 (x, ( 3 (t - h, x))(B t - B t - h ) 2 
+h 1+a R a (t - h, t, x), 
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for all < t — h < t < T, x £ M. d , which also holds on the set Q,'. Consequently, we obtain that 

G(t,x) -G(t-h,x) = L(x,( 2 (t-h,x))(B t -B t _ h ) + h 2a R a (t-h,t,x) 

= (D z G)(x,(2(t-h,x))G 2 (x,( 3 (t-h,x))(B t -B t _ h ) (5.6) 
+h 2a R a (t - h,t,x). 

In particular, by virtue of the Holder continuity of the Brownian motion (|4.1ip and the assumption 
(T3) we see from (15. 6p that 

G(x, ( 2 (t, x)) - G(x, C 2 (t - h, x)) = h a R a (t - h, t, x). (5.7) 

Similarly, we can also derive the following (recall (|5.5|) ) : 

| F(x,(2(t,x))-F(x,(2(t-h,x)) = h a R a (t-h,t,x); 
\ L(x,( 2 {t,x))-L{x,( 2 (t-h,x)) = h a R a (t-h,t,x), 

for all < t — h < t < T, \x\ < m, on the set il'. Now combining (|5.5p — (|5.8p . we obtain that, 
possibly on a different set fi, with P(£l) = 1, 

G(x,(2(t-h,x))(B t -B t „ h ) 
= G(x, C 2 (t, x)){B t - B t - h ) - (D z G)(x, ( 2 (t - h, x))G 2 (x, £ 3 (t - h, x))[B t - B t _ h f 
+h 3a R a (t - h,t,x) 

= G(x,Q 2 (t,x))(B t -B t _ h ) - (D z G)(xX2(t,x))G 2 (x,(3(t,x))(B t - B t ^ h f 
+h 3a R a (t - h,t,x). 

Moreover, rewriting ()5.8[) as F(x, ( 2 (t — h, x))h = F(x, ( 2 (t, x))h + h 1+a R a (t — h, t, x), and noting 
from (|5.8p (recall definition (|5.5p ) that 

(D z G)(x, ( 2 (t ~ h,x))G 2 (x, C 3 (t - h,x))(B t - B t _ h f 
= (D z G)(x,( 2 (t,x))G 2 (x,( 3 (t,x))(B t - B t _ h ) 2 + h ia R a (t - h,t,x), 

and that 

(D z G)(x, ( 2 {t - h, x))G 2 (x, ( 3 (t - h, x))h 
= (D z G)(x, C 2 (t, x))G 2 (x, C 3 (t, x))h + h 1+a R a (t - h, t, x), 

we obtain from (15.31) that 



((t,x)-({t-h,x) 

[F(x, C 2 (t, x)) - \(D z G){x, C 2 (t, x))G 2 {x, ( 3 (t, x))}h + G(x, C 2 (i, x))(^ - B t _ h ) 
-i(£> z G)(x, C 2 (t, x))G 2 (x, Ca(*, x))^ - B t „ h ) 2 + h 3a R a (t - h, t, x). 
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Finally, since a G (|, ^) is arbitrary, the proposition follows. ■ 
5.2 Backward temporal-spatial expansion 

We now give the complete statement of the backward temporal-spatial expansion. 

Proposition 5.2 Assume that ( G C°'( 3 )(F B , [0,T],M d ) satisfying (Tl)-(TS). Then, for all m G 
N and a G (|, ^) t/iere exists some subset $7 C S7 ,stic/i i/wi -P(S7) = 1, and i/iai on $7, /or all 
< t - h< t < T, x, k G R d , 

((t-h,x + k)-((t,x) = ah + b(B t „ h -B t ) + ^(B t _ h -B t ) 2 + (p,k)+^(Xk,k) 

+ ( q, k){B t _ h - B t ) + {h + \k\ 2 f a R a {t -h,t,x,x + k), (5.9) 

where 



a = 


- {F(x,( 2 (t,x)) -±(D z G)(x,(2(t,x))G 2 (x,(s(t,x))} 


b = 


G(x,<Z 2 (t,x)); 


c = 


(D z G)(x,( 2 (t,x))G 2 (x,( 3 (t,x)); 


P = 


(D x Q(t,x); 


X = 


D 2 x ((t,x); 


Q. = 


{D x G){x,( 2 {t,x)) + (D z G){x,( 2 (t,x))D x ( 2 (t,x). 



Proof. As in the forward expansion case, we need only show that for all m G N and a G (|, ^) 
there exists some subset $7 C $7 of full probability such that on 17, for all < t — h < t < T, 

x,k£ R d , 

C(t - h,x + k) - ((t - h,x) 
= D x ((t,x)k + ±{(D 2 x ((t,x)k,k) (5.11) 
+{(D x G)(x, ( 2 (t, x)) + (D z G)(x, ( 2 (t, x))D x ( 2 (t, x)}k(B t _ h - B t ) 
+(h + \k\ 2 f a R a (t -h,t,x,x + k). 

From the usual Taylor expansion with the remainder in the Lagrange form, we have 

C(t - h, x + k) - ((t - h, x) 
= D x ((t - h,x)k + -((DlC(t - h,x)k,k) + \k\ 3 R a (t -h,t,x,x + k), 
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for < t — h < t < T, x,k £ P-a.s. Moreover, from Proposition 15. II it follows that 



D x ({t,x)-D x ((t-h,x) 

D x {G{x, C 2 (i, x))]{B t - Bt- h ) + h 2a R a (t - h, t, x) 

{(D x G)(x, C 2 (t, x)) + (D z G)(x, ( 2 (t, x))D x ( 2 (t, x)}(B t - B t -h) + h 2a R a (t - h, t, x) 



and 



Dl((t, x) - Dl((t -h,x) = h a R a (t - h, t, x). 



Consequently, 



C(t - h,x + k) - C(t - h,x) 



D x ((t,x)k + ^((Dl((t,x)k,k) 

+{(D x G)(x, &(t, x)) + (D z G)(x, C 2 (t, x))D x ( 2 (t, x)}k(B t _ h - B t ) 

+ \k\ 6a R a (t - h,t - h,x,k) + h 2a kR a (t - h, t, x) + h a (R a (t - h, t, x)k, k). 



Thus, by virtue of 



kf a + h 2a \k\ + h a \k\ 2 < C m {h + \k\ 2 ) 



for h E [0,T] and \k\ < m, we derive (|5.1ip . This, combined with (|5.3p . leads us to the backward 



Finally, by combining Proposition 14.21 and Proposition 15.21 we have completed the proof of 
Theorem 12.31 

6 Application to Stochastic PDEs 

Having tried so hard to develop the various forms of stochastic Taylor expansion, as an application 
in this section we shall try to use it to study the stochastic viscosity solutions for fully nonlinear 
SPDE, following the idea that we developed in our earlier work [3]. In order not to over-complicate 
the computation we shall consider the following simpler version of the fully nonlinear SPDE (|2. 1[) : 



u(t, x) = u (x) + / f(x,u,Du,D 2 u)ds+ g(x, Du) o dB s , (t, x) G [0, oo) x M. d . (6.1) 



Compared to SPDE (I2.1j) . as well as Example 12.51 we see that here the diffusion coefficient g 
in (|6.ip is independent of u. The general case could be treated in a similar way but with more 
complicated expressions. Since our purpose here is to outline our idea of a new definition of 
stochastic viscosity solution, without adding too much technical complexity into this already 



temporal-spatial expansion (|5.9|) . 
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lengthy paper, we shall leave the study of the general case to a forthcoming paper. We should 
also note that in (|6.ip we are using the Stratonovich integral instead of the Ito integral for the 
simplicity of the presentation, the following relation is worth noting: 

g(x, Du) o dB t = g(x, Du)dB t + -D z g(x, Du)D x [g(s, Du)]dt. (6.2) 

It is worth pointing out that even in this simplified form, the nonlinearity of the function g on 
Du already makes it difficult to apply the rough path approach of [7] directly here. 

To explain our idea of the definition of stochastic viscosity solution, let us first apply Theorem 
12.31 to the regular solution u. Bearing the relation (|6.2p in mind we have: 

c 1 

u(t + h,x + k)-u(t,x) = ah + b(B t+h -B t ) + -(B t+h -B t ) 2 + { P ,k) + -(Xk,k) 

+ (q, k)(B t+h - B t ) + {\h\ + \k\ 2 f a R a>m (t, t + h,x, k), (6.3) 

where 

a = f(x,(u,Du,Du 2 )(t,x)), 

b = g(x,Du(t,x)) (6.4) 

c = (D z g(x,Du(t,x)),D x [g(x,Du(t,x))]), 

p = Du(t,x), X = D 2 u(t,x), 

Q = D x [g((x,Du)(t,x)\ = D x g(x,Du) + D 2 u(t,x)D z g(x, Du(t,x)). 

It is interesting to note that in this simple case the terms involving D 2 u can be written collectively 
as (suppressing variables): 

/(•••, D 2 u) + ([D 2 u]D z g(x, Du)(B t+h - B t ),k) 

+ l -{[D 2 u]D z g{x, Du),D z g(x, Du))){B t+h - B t ) 2 + ^({D 2 u)k, k) 

= /(•••, D 2 u) + - {[D 2 u](D z g(x, Du)(B t+h - B t ) + k),D z g(x, Du){B t+h - B t ) + k) . 

Compared to the classical deterministic Taylor expansion, and the SPDEs studied in [4] in 
which g is independent of Du, we can see that in a general case the terms involving D 2 u becomes 
much more complicated, and our previous method (via Doss-Sussmann) will face a fundamental 
challenge, especially in the uniqueness proof. We therefore will try to find a different approach to 
define the stochastic viscosity solutions, using the stochastic characteristics introduced by Kunita 
|12j . combined with our results on stochastic Taylor expansions. This new method also reflects the 
basic ideas of the works of Lions and Souganidis [T3] , [H] , [TB] , and [T7] , and in a sense includes 
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our previous work [3] as special case. For simplicity we shall now assume all processes involved 
are real valued. 

To motivate our definition of the stochastic viscosity solution let us suppose in a first step 
that the coefficients of SPDE (|1.4|) are sufficiently smooth and that this equation admits a regular 
solution u G C 0,oo (F s , [0, T] x R). Under this assumption, we compare the solution u with a 
smooth test field ip G C°'°°(F B , [0,T] x R) defined as the unique solution of the equation 

d(p(t,x) = 9(t, x)dt + g{x, Dcp(t, x)) o dB t , (6.5) 
cp(0,x) = <po(x), 

for t G [0,T], x G R d , where 9 G Cffi[0,T\ xR),g£ C^(R x R; R), and 990 G C~(R). 

Let us now fix an J-" B -measurable [0, T]-valued random variable r and an J 75 -measurable Re- 
valued random variable £. We say that u — ip achieves a local left-maximum in (r, £) if for almost 
all cj G {r < T} there is some p > (which may depend on u>) such that 

(u - <p)(u,t,x) < (it - ip)(u,T(u),^(u;)), 

for all (t,x) G [0,T] x R d with i G ((r(w) - p)+,r(o;)] and \x - < p. 

Following the approach by Kunita (see Theorem 6.1.2 in [12]), we formally introduce the 
following stochastic characteristics. For the moment let us assume that the random time r is 
actually deterministic to avoid further complications (keep in mind, however, that the Taylor 
expansion will hold even for the arbitrary random time r!). 

<j>t{x,z) = x- D z g((4> s ,Xs)(x,z)) o dB s , 



rj t (x,y,z) = y + J {g(((p s , Xs)(x, z)) - Xs(x, z)D z g(((p s , Xs)(x,z))} o dB s , (6.6) 
Xt(x,z) = z + J D x g((cf) s ,Xs)(x,z)) odB s , 
where t G [0, r], (x, y, z) G R 3 , and we hope to be able to define a transformation tp(t, x) by 

tp(t, 4> t (x, Dt/)(t, x))) = rft(x, if)(t, x),Dip(t, x)). (6.7) 
For notational simplicity let us now denote: 

l(x,z) = -D z g(x,z); h(x, z) = g(x, z) - zD z g(x, z)\ k(x, z) = D x g(x, z). (6.8) 
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Remark 6.1 We note that if the function g is linear in Du, then the situation will become much 
simpler, and the arguments below would become straightforward. To be more precise, let us 
consider the following two cases: 

(i) g{x,Du) = (H,Du), where H is a constant vector. Since D x g = D x V(x) and D z g = H, 
the Taylor expansion (|6.4p is drastically simplified. Also, the characteristics (|6.6p becomes almost 
triviel: x = z i V = Vi an d <fr = x — H{t — r). 

(ii) g(x, Du) = V(x)Du (see [7])- In this case one has D x g(x, z) = D x V{x)z, D z g(x, z) = V(x), 
thus the Taylor expansion (|6.4|) will also become much simpler. Furthermore, the characteristics 
(|6.6p now become "disentangled" SDEs: 

(pt(x) = x - J V((j) s {x))odB s , 
f]t(x,y,z) = y, 
Xt(x,z) = z+ D x V(4> s (x))x s {x,z) o dB s , 



and (|6.7p takes the special form: ip(t,x) = ip{t, (f>t(x)). In other words, the transformation tp is 
globally well-defined by an easy and explicit expression. ■ 

In what follows we shall denote O = (l,h,k). Also, for any function 7 = 7(x,y,z) we denote 
V7 = (-D27, -D y 7, D z r y) T . Then we can write, for example, the first equation in (|6.6p in the ltd 
integral from: 

4>t(x,z) = x+ l((<f> s ,x.s)(x,z)) o dB s (6.9) 



.>■+ I l(Vl,e)((^ s ,Xs)(x,z))ds+ I l{(4> s ,Xs){x,z))dB s 



2 

Now, treating the random fields r/ and x i n the same way, and applying the stochastic 

3' 2/' 



(backward temporal) Taylor expansion to (4>, r/, x) we have, for all a E ( 1 1 



<k(x, z) = x + l(x, z)(B t - B T ) + \ ( V/, 9 )(x, z)(B t - B T f + \t - r\ 3a R at7n , 
r, t (x, y,z)=y + h{x, z)(B t - B T ) + - ( Vh, G )(x, z){B t - B T ) 2 + \t- T\ 3a R a , m , (6.10) 
Xt (x, z) = z + k{x, z){B t -B T ) + -( Vfc, 6 >(x, z){B t - B T f + \t- T\ 3a R a , m . 
On the other hand, writing (|6,5p in an ltd integral form we have 

f* 1 

ip(t,x) = ipo(x) + J {0(s,x) + -D z g(x,Dip(s,x))D x [g(x,Dtp(s,x))}}ds 

+ / g(x,D<p(s,x))dB s (6.11) 



A 



fo(x) + J F(s,x)ds + J G(x,(2(s,x))dB s 
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where F and G are denned in an obvious way, and with 

(i(t,x) = <p(t,x), ( 2 (t,x) = D<p(t,x), (3(t,x) = (D(p,D 2 <p)(t,x). 

Note that differentiating the both sides of (|6.1ip we have 

dD<p(t, x) = D x F(t, x)dt + {D x g(x, Dcp(t, x)) + D z g(x, D<p(t, x))D 2 x Lp(t, x))}dB t . (6.12) 

Applying the backward temporal Taylor expansion (Theorem 2.3) again to the random field (p 
around any point (r, <j), with 

F 2 (t,x) = D x F(t,x); G 2 (x,z 1 ,z 2 ) = D x g{x,zi) + D z g(x,z 1 )z 2 , 

we obtain, after some simple cancelations, that 

<p(t,x) = ip{T,q)+e{T^){t-T)+g{s : Dip{r,q)){B t -B T ) 

+^{{D z9 D x g}^, D^{t, 0) + {D z g{q, D<p(r, ?))) 2 Z?^(r, s)}(B t - B T f 
+D<p(r, q)(x - ?) + D x [g(; D<p(; -))](t, s)(x - s)(B t - B T ) (6.13) 

+ ~Dlv(T, S)(X - ?) 2 + (\t -T\ + \X- s\ 2 f a Ra,m. 

We note that the above holds for all (r, <;) and all oj € $l a ,m- 

Let us now define the desired random field if). We should note that our main purpose here is 
to find such a transformation so as to eliminate the stochastic integral. In other words, we shall 
look for such if) that has the following first order Taylor expansion: 

f lf){t,x) =<p{T,x)+d 1 Tllj{T,x)(t-T) + \t-T\ 3a R c 



„ ., (6-14) 

{ Dtp(t,x) = D<p(T,x) + d^Dip(r,x)(t - r) + \t - r\ ia R a ^ m . 

Here 9 t ~ denotes the left partial derivative with respect to t. To this end, we note that (|6.6p and 
(|6.7p imply that for t = r, one has 

ip(r, x) = tp(r, 4> T {x, Dtp(r, x))) = r) T (x, ip(r, x),Dip(r, x)) = ip(r, x), x G R, (6.15) 

and hence Dip(r, x) = Di/j(t, x) holds for all x 6 R as well. Next, we look at the Taylor expansion 
for both cf) and r]. Recalling (j6.8|) and the Taylor expansions (|6.10p . We shall first take x = 
(x,ifi(t,x),Dip(t,x))), and then replacing x by (f>t(x, Dif)(t,x)) and ? by x in (|6.13p . It should be 
noted that after these substitutions the remainder will look like R atm (t,r,x,4>t(x,Dif>(t,x))), a 
slightly more complicated form than the original ones. To make sure the accuracy of the expansion 
in what follows we shall denote 

Ra,m,loc= SUp _ \R a ,m(t, S, X, <j) t (x, Dlf)(t, x))) \ . 

t,se[0,T]; xeBm(O) 
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Then it is not hard to show that there exists an increasing sequence {&£,£ > 1} C J~ B with 
lim^+oo P(Qe) = 1, such that, 

Keeping such a modification in mind, we can now proceed to write down the Taylor expansion: 

(p(t,(j)t(x,Di()(t,x))) 
= ip(r, x) + 0{t, x)(t - t) + g(x, Dy(r, x)){B t - B T ) 

+\{[D z gD x g]{x, D<p{r, x)) + (D z g(x, D<p(r, x)) 2 D 2 x ip( T , x)](B t - B T ) 2 
-D<p(r, x)D z g(x, D<p(r, x)){B t - B T ) 
-D<p(T,x)D 2 z g(x,Dtp(T,x))dt Di>(r,x)(t - r)(B t - B T ) 

+^£V(t, x) ( V/, 8 )(x, D<p(r, x)){B t - B T ) 2 (6.16) 
-^D 2 Mr,x)(D z g(x,D^(T,x))) 2 (B t - B T ) 2 

-D x g(x, D<p(r, x))D z g(x, D<p{r, x))(B t - B T ) 2 + \t- T\ 3a R a ^ loc 
= ip(r, x) + 6{t, x)(t - r) + g{x, D(p(r, x)){B t - B T ) 

-^[D z gD x g}(x,Dip(r,x))(B t - B T ) 2 - D<p(t,x)D z9 (x, D<p(T,x))(B t - B T ) 

+^D^(r, x) ( V/, )(x, Dcp(T, x))(B t - B T ) 2 

-Dtp(T,x)D 2 g(x 1 Dip(T,x))di DiJj(T,x)(t - r)(B t - B T ) + \t - r\ 3a R a , m , loc . 
On the other hand, from (|6,6p we see that 



fjf(t,x) = ri t (x,ip(t,x),D4>(t,x))-h(x,Dip(t,x))(B t -B T ) 

~{Vh,@) (x, DrP(t, x))(B t - B T ) 2 + \t - r\ 3a R a , m> i oc . 

Using the form (|6.14p and the smoothness assumptions on all the coefficients, and then noting 
(|6,15p one can easily rewrite above as 



ip(t,x) = rj t (x,t/j(t,x),Difj(t,x)) - h(x,Dip(T,x))(B t - B T ) 
-D z h(x,D<p(T,x))diDi;(T,x)(t - r)(B t - B T ) 
-i ( Vh, Q)(x, Dip(r, x))(B t - B T ) 2 + \t- r\ 3a R a ^ loc . 
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Now recalling the relation (|6,7p . the fact (|6.16p . and the definition of h (|6.8p . we obtain that 

tp(t,x) = T]t(x, ip(t, x), Dip(t, x)) — h(x, D(p(r, x))(Bt — B T ) 

~ ( Vh, 9 )(x, D<p{r, x))(B t - B T ) 2 + \t- r\ 3a R a:m 
= <p(r, x) + 0(t, x)(t - r) + g(x, Dip{r, x)){B t - B T ) 

-^[D z gD x g]{x,D^{T,x)){B t - B T f - D V (T,x)D z g(x, D(p(r,x))(B t - B T ) 

+^D^{t, x) ( VI, )(x, D<p(r, x)){B t - B T f 

-Dip(r, x)D 2 x g{x, D<p(r, x))9fDip(r, x)(t - r)(B t - B T ) (6.17) 

-[g{x, Dtp(r, x)) - Dip(r, x)D z g(x, Dtp(r, x))](B t - B T ) 

-D z h(x, D<p(T, x))d^D^(r, x)(t - T){B t - B T ) 

-i ( Vh, Q)(x, D<p(r, x))(B t - B T f + \t- r\ 3a R atm>loc 

= ip(r, x) + 6(t, x)(t - r) - ^[D z gD x g](x, D<p(r, x)){B t - B T ) 2 

+^D<p(r, x) ( VI, 8 )(x, D<p(t, x))(B t - B T f 

-D<p(r, x)Dlg(x, D<p(r, x^D^t, x)(t - r){B t - B T ) 

-D z h(x, D<p(T, x))d^Di,{r, x)(t - r)(B t - B T ) 

- X - ( Vh, 6 ){x, Dip{r, x))(B t - B T f + \t - r\ Za R a ^ loc . 

Since V/ = (—D xz g,0, —D 2 g) and Vh = (D x g — zD xz g,0, —zD 2 g), one can check that 

~{D z gD x g] (x, D<p(r, x)) + ^D<p(r, x) (VI, Q) (x, D<p(r, x))-^{Vh,Q) (x, D<p(r, x))=0, 

and 

-Dipir, x)B\g(x, Dip(r, x))6£ Zty(r, x))(t - r)(B t - B T ) 

-D z h(x,D(p(T,x))dfDi/>(r,x)(t - r){B t - B T ) = 0. 

In other words, f|6. lTj) leads to that 

ij){t, x) = ip{r, x) + 9{t, x)(t - T ) + \t- T\ Za R a ^i oc . (6.18) 

Clearly, (|6. 18[) indicates that 

dti){T,x) = 9{t,x), dtDip{T,x) = D9{t,x). 

We can now define the notion of stochastic viscosity solution for (|1.4p . 
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Definition 6.2 A random field u € C([0,T] x M. d ) is called a stochastic viscosity subsolution 
(resp. supers olution) of if u(0,x) < uq(x) (u(0,x) > uq(x)) for all x £ M. d and if for 

any r 6 L°{F B , [0, T]), £ G L°(T B ,R d ), and any (not necessarily adapted) random field if G 
C°' 2 ([0,T] x M d ) having the expansion \6.13\) . it holds that 

e(r,0 <(>) f(t,u(T,0,Dtp(T,0,D 2 MT,0). 

P-a.s. on the subset of Vt on which u — <p < (>)(u — <p)(t,£) = at a left neighborhood of (r, £). 
A random field u 6 C([0,T] x M. d ) is called a stochastic viscosity solution of \l-4\) > if it is both a 
stochastic viscosity subsolution and a supersolution. 

Remark 6.3 We would like to note that in Definition 16.21 the (possibly anticipating) test func- 
tions (p can be directly defined by the expansion (|6.13p with <p(t,£) being replaced by u(t, £) 
and (9(t, £), D x ip(r, £), D 2 (p(r, ^)) being replaced by a triplet of random variables (f3,p,A). The 
main advantage here is that the random field ip is now defined globally, overcoming the essential 
difficulties in the theory of stochastic viscosity solution thus far, and will significantly facilitate 
the uniqueness proof. Moreover, on the subset of Vt where u — <p achieves the local left-maximum 
(resp., left-minimum) at the point (r, £), the triplet (/3,p, A) can be considered as a stochastic 
sub- (resp. super- )jet, as it was traditionally done. These issues will be further explored in our 
forthcoming publications. ■ 

In the rest of this section we shall verify that a regular solution must be a stochastic viscosity 
solution in the sense of Definition 16.21 which will provide a justification for our new definition. To 
this end, let us assume that the coefficient / is proper in the following sense 

(HI) The function F(t,u,p, X) = —f(x,u,p,X) is "degenerate elliptic". That is, f is con- 
tinuous in all variables, and is non- decreasing in the variable X . 

Assume that u is a regular solution to (jl.4p . then we have 

u(t,x) = u (x) + [ F(s,x)ds+[ g(x,Du(s,x))dB s , t>0, (6.19) 
Jo Jo 

where 

F{t, x) = f{x, (u, Du, D 2 u)it, x)) + -D z g(x, Du(t, x))D x [g{; Du(t, •))] (t, x). 

For any given pair of random variables (r, £) and an arbitrary test field (p such that u — <p attains 
a local left-maximum at (r, £) on a subset of Q with positive probability, let ip be the process 
associated to (p by (|6.18p (where the remainder can be chosen to be zero!), and <p be defined by 
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We first apply the Taylor expansion on u at point (r,x), and evaluated at (t, (j>t(x, Dij)(t, x))) 
and then use the expansion of (f>t(t,Dip(t,x)) (recall (|6.1U|) ) to get 

u(t,(p t (x,Dip(t,x))) 
= u{t, x) + f(x, (u, Du, D 2 u)(t, x))(t - r) + g(x, Du(t, x)){B t - B T ) 

+- {D z g(x, Du(t, x))D x g(x, Du(t, x)) + D z g(x, Du(t, x)) 2 D 2 xx u(t, x)} (B t - B T f 

+D x u(t, x)(4> t (x, Dip(t, x)) - x) + ^D 2 xx u(t, x)(4> t {x, Dtf;(t, x)) - x) 2 
+D x [g(x, Du(t, x))]{B t - B T )(<f> t (x, Lty(t, x)) - x) + \t - T\ 3a R^ m (r, t, x, (f> t (x, Di/)(t, x))) 
= u(t, x) + f(x, (it, Du, D 2 u)(t, x)(t — r) + g(x, Du(t, x)){B t — B T ) 

+^{[D z gD x g](x,Du(T,x)) + (D z g(x,Du(T,x)) 2 D 2 x u(T,x)}(B t - B T ) 2 (6.20) 

-Du(t, x)D z g(x, D<p(r, x)){B t - B T ) 

-Du(t, x)D 2 g(x, Dip(r, x))dfDij>(T, x))(t - r)(B t - B T ) 

+^Du(r, x) (V/, 0) (x, D<p(r, x))(B t - B T ) 2 

+^D 2 x u(t, x)(D z g(x, D<p(r, x))) 2 {B t - B T ) 2 
-(D x g)(x, Du(r, x))D z g(x, D<p(r, x))(B t - B T ) 2 

-(D z g)(x, Du(t, x))D z g{x, Dip{r, x))D 2 x u{t, x)(B t - B T ) 2 + \t- r\ 3a R a>m>loc . 

We will now replace x by the random point £. To specify the "left-neighhood" required for the 
viscosity property, we also define for any p > the following subset of $7: 

r r,'| = {w : (n - (p)(oj,t,x) < (u - (p)(lj,t(<jj),€(<jj)), (t - t(u))~ < p and x G B p (£(u>))}, 

To wit, T^'j? is the subset of on which u — ip attains a left local maximum at (r, £). Now setting 
x = £ in ()6.20p . and noting that Du(t,£) = D(p(r,£) and u(r, £) = ^(t, £), we have 

u(t,MZ,Dij(t,0)) 
= <T, + f(£, (u, Du, D 2 u)(t, - r) + g(£, Du(r, 0)(B t - B T ) 

-~[D z gD x g](t;,Du(T,0)(Bt - B T ) 2 - Du(r,x)D z g((, Du{r,i)){B t - B T ) 

+^Du(t, e)(V/, 0>(z, Du{t, 0)(B t - B T ) 2 

-Du(T,OD 2 g(Z,Du(T,0)diD^(T,0)(t - r)(B t - B T ) + \t- r\ 3a R a ^ loc , 
almost surely on . Consequently, from (|6.16p we obtain that, P-a.s. on T^, 

o > {u(t,Mt,D^(t,m-<T,o) -&(t,Mt,D^(t,m-v(T,o) 

= {/(£, (u, Du, D 2 u)(t, 0) ~ 0(t, 0}(t - t) + \t - r\ 3a R a ^ loc . 
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Since K r on T^, we deduce that 

/(£, (u, Du, D 2 u)(r, 0) ~ 0(t, > 0, P-a.s. on 1^. 

Finally, since the mapping X i-> f(x,u,p,X) is non-decreasing, thanks to (HI), and since 
£>V(t, > D 2 u(t, on we have 

/(£, L^, D\){t, 0) > 0(r, e) P-a.s. on 1^. 

This proves that the classical solution u is a stochastic viscosity subsolution. That u is also 
a supersolution can be proved using a similar argument. Therefore, a is a stochastic viscosity 
solution. 
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